New Nuclear Physics for Big Bang Nucleosynthesis 



O 

(N 

< 

O 

u 

: 

Or 

i < 

o ■ 

c3 



> 

00 

oo 
q 

00 

o 
o 



X 



Richard N. Boyd 1 , Carl R. Brune 2 , George M. Fuller 3 , and Christel J. Smith 4 
1 Lawrence Livermore National Laboratory, Livermore, CA 94551 
2 Department of Physics and Astronomy, Ohio University, Athens, OH 45701 
3 Department of Physics, University of California, San Diego, La Jolla, CA 92093-0319 
^Physics Department, Arizona State University, Tempe, AZ 85287-1404 
(Dated: August 5, 2010) 

We discuss nuclear reactions which could play a role in Big Bang Nucleosynthesis (BBN). Most 
of these reactions involve lithium and beryllium isotopes and the rates for some of these have 
not previously been included in BBN calculations. Few of these reactions are well studied in the 
laboratory. We also discuss novel effects in these reactions, including thermal population of nuclear 
target states, resonant enhancement, and non-thermal neutron reaction products. We perform 
sensitivity studies which show that even given considerable nuclear physics uncertainties, most of 
these nuclear reactions have minimal leverage on the standard BBN abundance yields of 6 Li and 
7 Li. Although a few have the potential to alter the yields significantly, we argue that this is unlikely. 

PACS numbers: 14.60.Pq; 14.60.St; 26.35,+c; 95.30.-k 



I. INTRODUCTION 

The basic details of big bang nucleosynthesis (BBN) 
have been well understood for several decades [lpfl5| , and 
nucleosynthesis predictions are in reasonable agreement 
with the abundances of 2 H, 3 He, and 4 He observed in 
metal poor environments. However, there are new chal- 
lenges in achieving agreement between theory and ob- 
servation for the lithium isotopes, 6 Li and 7 Li. These 
challenges stem from recent observations of metal-poor 
halo stars [l6| and the high precision measurement of 
the baryon-to-photon ratio of the Universe by WMAP 
[l7l - l2l1 |. The WMAP value, which confirms earlier infer- 
ence of this quantity from deuterium [22], [H[ , is likely to 
be improved on by future cosmic microwave background 
observations (e.g., Planck [H|). 

The abundance of 7 Li has been a longstanding prob- 
lem, as the predicted abundance is higher than observed 
in metal poor halo stars 25] by roughly a factor of three 
at the WMAP baryon-to-photon ratio. Attempts have 
been made to reconcile this difference, primarily on the 
basis of stellar processing of 7 Li. Although it can be ar- 
gued that some 7 Li would be destroyed by the stars in 
which it is observed, it has been found to be very dif- 
ficult to justify enough destruction to bring theory and 
observation into agreement [26[ . 

Recently Very Large Telescope observations [l6| sug- 
gested an abundance of 6 Li on the surface of some metal- 
poor halo stars roughly ~ 1/30 of that of 7 Li. Although 
this observation has been challenged by other work [27j . 
if it is correct, it represents a 2-3 order of magnitude 
discrepancy between theory and observation. Much is 
riding on this discrepancy. If 6 Li must be made by post- 
BBN non-thermal processes, then the abundance of this 
nucleus conceivably may be an indirect probe of new 
physics , specifically heavy particle decay and dark matter 

Nollett et al. [3l[ have provided an insightful and 
comprehensive study of the nuclear reactions which bear 



on 6 Li production and destruction, while Cyburt and 
Pospelov [32| have explored the necessary issues for a 
pure nuclear physics solution to the 7 Li overproduction 
problem. In this paper we explore related issues, find no 
clear paths to solution of either lithium problem, but do 
provide rates for incorporation into BBN codes as well 
as studies of the sensitivity of lithium isotope production 
and destruction to uncertainties in these rates. It should 
be kept in mind that although these rates may make lit- 
tle difference in standard BBN abundance yields, using 
BBN to constrain non-standard scenarios with new par- 
ticle physics may require new levels of precision in some 
parts of the BBN nuclear reaction network [33|, HH • 

Although the standard hot BBN model code contains 
most of the reactions that could be relevant to BBN, 
those involving reactions on short-lived nuclei are not 
well studied, and in some cases have not been included. 
Given the current level of precision of BBN calculations 
[Tol 135 . 36] , it seems appropriate to reexamine the BBN 
network to be sure that all possible reactions are in- 
cluded, and to study the potential effects of those re- 
actions for which data do not exist. Candidate reactions 
here include 7 Be( 3 H, 4 He) 6 Li and, possibly, other reac- 
tions on 7 Be, particularly because mass-7 is made pri- 
marily as 7 Be at the WMAP baryonic density. Another 
reaction that might be of particular interest for synthe- 
sizing 6 Li could be 3 H( 3 He,7) 6 Li, which was measured, 
but in a very difficult experiment. 

Another aspect of BBN that we have studied is the in- 
clusion of nonthermal particles. It has generally been as- 
sumed that the BBN reactions occur in an environment 
where ions have thermal, Maxwellian velocity distribu- 
tions. This approximation could be violated at some 
level by reactions induced by highly-energetic particles 
produced by exothermic reactions. Again, in the interest 
of studying the precision that can be expected of BBN 
calculations [13, HI, Hfjj], we have studied the effects of 
non-thermal particles. The most likely scenarios for this 
process involve neutrons produced in 3 H(d, n) 4 He and 
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the protons produced in 3 He(c?,p) 4 He. These reactions 
involve relatively abundant nuclei in the entrance chan- 
nel, have large cross sections, and produce nuclei with 
energies in excess of 10 McV. Some aspects of this were 
studied recently 37[. That work showed that thermal- 



ization of the non-thermal charged particles occurs suffi- 
ciently rapidly, through electromagnetic processes, that 
those particles have only a tiny effect on the calculated 
BBN abundances. However, that study only investigated 
the thermalization of protons and did not consider neu- 
tron thermalization, which occurs over a much longer 
timescale than for protons. It might be thought that 
the non-thermal neutrons could affect the BBN abun- 
dances, specifically that of 7 Be, through 3 He(n,p) 3 H and 
through neutron-induced reactions on 7 Be. The former 
reaction could reduce the abundance of 3 He, and there- 
fore, of 7 Be, since it is made from 3 He( 4 He,7) 7 Be. The 
latter reactions might reduce the abundance of 7 Be after 
it was made, and 7 Be(n, 2 H) 6 Li, which is endothermic, 
could produce 6 Li. Therefore we have considered the ef- 
fects of non-thermal neutrons on the 7 Be abundance, and 
on possible breakup of 2 H as well. 

The lowest-lying first excited states of the lightest nu- 
clei are those of the mirror nuclei 7 Li and 7 Be. The exci- 
tation energies of these excited states are within a factor 
of 3 or 4 of relevant BBN thermal energies. Accordingly, 
we have examined how potentially key nuclear reactions 
like 7 Be(d,7) 9 B, 7 Be(d,p)2a, and 7 Be(d, 3 He) 6 Li could 
be altered if they were to proceed through a thermally 
populated first excited state, 7 Be*. Cyburt, Ref. [32|, ar- 
gued that resonant enhancement of 7 Be(<i, p)2a could in 
principle solve the first of the lithium problems outlined 
above, potentially reducing the 7 Li BBN yield by a fac- 
tor of 3 or so. We revisited that possibility, but also find 
that if this reaction were to proceed through 7 Be* a dif- 
ferent resonant enhancement channel can come into play. 
In the end, however, standard BBN yield alterations are 
found to be small due to the small thermal population of 
7 Be*. 

In what follows we discuss the BBN reaction network 
in section II and the results of our BBN calculations with 
the new rates in section III. In section IV we discuss non- 
thermal neutron effects and in section V we discuss ther- 
mal population of excited states and resonant enhance- 
ment. We give conclusions in section VI. 



have not generally been included. In most cases this is 
because these reactions are difficult to study, for exam- 
ple, because either the target nucleus or the beam nu- 
cleus, or in some cases both, are radioactive. The inputs 
to the BBN code, that is, the information from the ex- 
periments, is in the form of thermonuclear reaction rates, 
{av), defined as [381 ] 



(av) 



7T/i 



1/2 



fcT]- 3 / 2 J a(E)Ee- E / kT dE, 



(1) 



where \i is the reduced mass and E is the center of mass 
energy of the reactants. Nuclear astrophysical measure- 
ments usually produce an astrophysical S- factor, S(E), 
defined as 



S{E) = Ea(E)e 



bE 



-1/2 



(2) 



where e 



is (approximately) the Coulomb barrier 



tunneling probability, b — (2/i) 1 / 2 Tre 2 ziZi/h, and z\ and 
Z\ are the charge numbers of the entrance channel nu- 
clei. In the absence of resonances, S(E) will be fairly 
constant, and the thermonuclear reaction rate (av) times 
Avagadros constant, N A , can be shown [38[ to be: 



N A (av) = 4.34 x 10 5 



S(E) 
[fxziZi 



tV 



cm" 



mole ' 



(3) 



where r = 4.246[z 2 ZjV/T 9 ] 1 / 3 , with T 9 being the tem- 
perature in billions Kelvin, and S(E) is in units of keV- 
barns. Since most of the changes we made to the BBN 
reaction rates involved reactions that produced high ex- 
citation energies in the compound nucleus, we found that 
the S-factors we calculated were usually sufficiently non- 
resonant that we could use the above prescription. 



II. BBN NUCLEAR REACTION NETWORK 



In this paper, we discuss several reactions that should 
naturally be a part of the BBN reaction network, but 



In one case we needed to calculate the reaction rate for 
a narrow resonance, for which the formula is 38] 



J 



N A (av) = 1.535 x 10 12 (w7) 3 /x- 3/2 T 9 ~ 3/ V 1L605 ^ 6/T9 cm 3 s" 1 mole" 



(4) 



where (^7)3 = wril^/T (given in keV), with Ti and channels, T is the total width, and u> is the spin statistics 
T2 the resonance partial widths in the entrance and exit 
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Reaction Key 




n -> J H 


1 


x H(n,g) 2 H 


2 


2 H( 2 H,n) 3 He 


3 


z H( 1 H,g) 3 He 






4 


2 H(n,g) 3 H 




3 H( 2 H,n) 4 He 


5 






3 H( 4 He,g) 7 Li 


6 


He(n, H) H 


7 


3 He( 2 H, 1 H) 4 He 


8 


3 He(n,g) 4 He 




3 He( 4 He,g) 7 Be 


9 


7 Li( 1 H, 3 He) 4 He 


10 


7 Be(n, 1 H) 7 Li 


11 


7 Be(n, 2 H) 6 Li 


12 


7 Li(n,2n) 6 Li 
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4 He( 2 H s) 6 l\ 


14 


6 Li( 1 H, 4 He) 3 He 


15 


7 Li( 3 H,n) 9 Be 


16 


7 Li( 3 He, 1 H) 9 Be 




7 Li( 2 H,n) 8 Be 


17 


7 Li( 1 H,g) 8 Be 




7 Li( 3 He, 2 H) 8 Be 




8 Be -> 4 He+ 4 He 




7 Li( 3 H, 4 He) 6 He 


18 


6 He-> 6 Li or 4 He+ 2 H 


19 


7 Li( 3 He,n) 9 B 


20 


3 He( 3 H,g) 5 Li 


21 


7 Be( 1 H,g) 8 B 


22 


7 Be( 3 H, 1 H) 9 Be 


23 


7 Be( 2 H, 1 H) 8 Be 


24 


7 Be( 3 H, 2 H) 8 Be 




7 Be( 3 He, 1 H) 9 B 


25 


7 Be( 3 H,n) 9 B 




9 Be( 1 H, 4 He) 6 Li 


26 



FIG. 1: Reaction network for BBN, as modified from [jj^]. The nuclei indicated in dotted boxes are completely unstable, 
whereas those in dashed boxes have some stable states and some others that are relevant to BBN that undergo particle decay 
to some other nuclei. Red lines indicate reactions newly added to the BBN code, and blue lines indicate reactions studied in 
the context of non-thermal neutrons. Double thickness lines indicate more than one possible reaction. 



factor = (2J+l)/[(2Ji + l)(2J 2 + l)], with J the total res- 
onance spin and J% and J2 the spins of the two entrance 
channel particles. Eq is the center of mass resonance en- 
ergy in the entrance channel in MeV. This formula does 
not apply to cases in which there are broad resonance 
structures; then the reaction rates need to be calculated 
by numerical integration or approximated in other ways. 

The study by Serpico et al. [5(| investigated the sensi- 
tivity of BBN abundance yields to most of the usually 
included reactions. We have added several "deuteron 
transfer" reactions. These can be quite strong, and con- 
ceivably could have an appreciable effect on the 7 Li and 
7 Be abundances. For example, the reaction rate for 
7 Li( 3 H,n) 9 Be, measured by Brune et al. (39|, was found 
to proceed to both the 9 Be ground state and to excited 



states. Since the rate to the ground state is 0.084 times 
the rate to the continuum states, destruction of 7 Li is 
far more likely than production of 9 Be by this reaction. 
We have also added to the code several reactions that 
occur on 7 Be, e.g., 7 Be(p, 7 ) 8 B H|, 7 Be(d,p) 8 Be and 
7 Be( 3 H,n) 9 B 42j, all of which create nuclides that un- 
dergo particle decay, and so result in destruction of 7 Be. 

Another reaction that has interesting potential is 
3 He[ 3 H,7) 6 Li, which has been studied experimentally 
[471 l48j |. has a huge positive Q-value (15.795 MeV), and 
can go to many excited states of 6 Li. The 6 Li ground 
state is stable, the second excited state decays to the 
ground state, and all other states undergo breakup to 
4 He+ 2 H. Thus this reaction can produce some 6 Li. Its 
effects were studied theoretically by Fukugita and Kajino 
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TABLE I: 
Reactions Added to the BBN Code 



Reaction 


Q-value (MeV) 


Comment 


Effect 


7 Li( 3 H,n) 9 Be(ground state) 


10.439 


1 (Ref. [39]) 


none 


7 Li( 3 H,n) 9 Be(excited states) 


10.439 


1 (Ref. [3911 


none 


7 Li( 3 He,p) 9 Be(ground state) 


11.202 


2 (Ref. [40, 411) 


none 


7 Li( 3 He,p) 9 Be(excited states) 


11.202 


2 (Ref [40 411) 


none 


7 Li( 3 He,n) 9 B 


9.352 


9 (Ref. [42]) 


none 


7 Li( 3 He, 2 H) 8 Be 


17.608 


8 (Ref. [43. 441) 


none 


7 Li( 2 H,n) 8 Be 


15.031 


8 (Ref. [43. 441) 


none 


7 Li( 3 He, 4 He) 6 Li(ground state) 


13.328 


4 (Ref. [451) 


none 


7 Li( 3 He, 4 He) 6 Li(excited states) 


13.328 


4 (Ref. [451) 


none 


7 Li( 3 H, 4 He) 6 He 


9.838 


3 (Ref. [461) 


none 


7 Be( 3 H, 4 He) 6 Li(ground state) 


14.208 


4 


small 


7 Be( 3 H, 4 He) 6 Li(excited states) 


14.208 


4 


small 


7 Be( 3 H,p) 9 Be(ground state) 


12.082 


6 


small 


7 Be( 3 H,p) 9 Be(excited states) 


12.082 


6 


small 


3 He( 3 H,7) 6 Li(low lying states) 


15.795 


5 (Ref. [47. 481) 


large 


9 Be(p, o;)^Li(grouiid state) 


2.126 


Ref. [49J 


none 


7 Be(p, -7) S B 


0.137 


Ref. [49J 


none 


7 Be( 2 H,p) 8 Be 


16.674 


8 (Ref. [43. 441) 


large 


7 Be( 3 H, 4 He) 6 Li 


14.208 


4 


small 


7 Be( 3 H,n) 9 B 


10.232 


7 


small 


7 Be( 3 H, 2 H) 8 Be 


12.641 


8 (Ref. [43,44J) 


small 


7 Be( 3 He,p) 9 B 


10.995 


7 


none 


7 Be( 2 H, 3 He) 6 Li 


-0.112 


10 


none 


7 Be*( 2 H, 3 He) 6 Li 


0.317 


10 


large 


7 Be*( 2 H,p) 8 Be 


17.103 


8 (Ref. [43,44]) 


large 


7 Be*(p, 7 ) 8 B 


0.566 


Ref. [49J 


large 



[5lj|, an d by Madsen [^3, and were found to contribute 
little to 6 Li production. However, there is some question 
as to the normalization of the cross section for that re- 
action; this will be discussed below. Another promising 
reaction is 7 Li( 3 He, 4 He) 6 Li [45[, which proceeds through 
high- lying resonances in 10 B, and may be able to pro- 
duce some 6 Li. Finally, 7 Be( 3 H, 4 He) 6 Li can both de- 
stroy mass-7 nuclei and produce 6 Li. This reaction has 
not been studied experimentally, since both nuclei in the 
entrance channel are radioactive. It has a large Q-value, 
so high-lying states in 6 Li that it can populate must be 
considered in the context of mass-7 destruction. 



We have terminated our reaction network at mass 9. 
We believe this is justified by the 9 Be(f>, a) 6 Li reaction 
rate [49(; its magnitude insures that nearly all of the 
9 Be made by 7 Li( 3 H,n) 9 Be or any other reaction will be 
returned to 6 Li. As pointed out by Boyd and Kajino [53j j . 
9 Be( 3 H,n) n B also has the potential to destroy 9 Be, but 
its rate is so much smaller than that for 9 Be(p, a) 6 Li [49] 
that it is of little consequence for destroying 9 Be, and 
results in so little production of mass- 11 nuclei that they 
will be far below observational limits. 



In the end, however, our attempts to augment the pro- 
duction of 6 Li were undone, principally by 6 Li(p, a) 3 He, 
the large cross section for which results in the destruction 
of nearly all of the 6 Li made in BBN. This was nicely set 
out in Nollett et al. [9]. In the next section we will add 
to this discussion. 



A. Comments on the Added Nuclear Reaction 
Rates 

Some comments are in order concerning the reaction 
rates we used in our calculations. These reactions are 
listed in Table I. 



1. 7 Li( 3 H,n) 9 Be (ground state) and 7 Li( 3 H,n) 9 Be(excited 
states) 

These rates have been determined experimentally [39l | , 
so are known accurately. The reaction to the ground 
state gives one neutron, whereas those to all the excited 
states give two neutrons due to the breakup of all 9 Be 
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excited states into two a-particles and a neutron. The 
relevant rates 39] are: 

N A (av{ 7 U( 3 H, n) 9 Be(ground state))) = (5) 
2.98 x 10 10 r 9 - 2/3 e- 11 - 327T 9" I/a 
x[l - 0.122T 2/3 + 1.32(T 9 4/3 0.127T 9 2/3 + 0.0742)- 1 ] 
cm 3 s _1 mole -1 , 

A A (crt;( 7 Li( 3 H,n) 9 Be(excited states))} = (6) 
11.9 x iV j4 (crw( 7 Li( 3 H,n) 9 Be(ground state))). 

2. 7 Li( 3 He,p) 9 Be (ground state) and 7 Li( 3 He,p) 9 Be (excited 

states) 

The first reaction was studied at higher energies [ioj . 
and extended to lower energies [4l[ , and an S-factor was 
determined for the reaction to the 9 Be(ground state). Al- 
though the S-factor exhibits a broad resonance structure, 
the resonances exist at relatively high energies. At en- 
ergies below 100 keV, where this reaction is relevant to 
BBN, the S-factor is essentially flat at 5.3 MeV-barn. In 
addition, from the spectrum [40j it can be inferred that 
the reaction to excited and continuum states is roughly a 
factor of five greater than that to ground state. Therefore 
the assumed rates, in units of cm 3 s _1 mole -1 , are: 

iV A (cru( 7 Li( 3 He,p) 9 Be(ground state))) = (7) 
5.87 x i iOT 9 - 2/3 e- 17 - 980T 9 " 1/3 

iV j4 (crw( 7 Li( 3 He,p) 9 Be(excited states))) = (8) 
2.94 x i iiT 9 - 2/ V 17 - 980T o~ 1/3 

3. 7 Li( 3 H, 4 He) 6 He 

Although this reaction makes 6 He, its ground state /3- 
decays with a half-life of 807 ms to 6 Li. It was studied 
at low energies [III, both to the ground state, and to 
the first excited state, which breaks up into 4 He and two 
neutrons. The S-factor to the 6 He ground state was esti- 
mated from data in Ref. [46[ to be 0.84 MeV-barn, and 
that to the first excited state to be 13.9 MeV-barn. For 
the two rates associated with this reaction, this gives 

V A (crt>( 7 Li( 3 H, 4 He) 6 He(ground state))} = (9) 
7.4 x lQ 9 T 9 - 2/3 e- u - 327T 9" 1/a cm 3 s" 1 mole" 1 , 

N A {av( 7 U( 3 U, 4 He + 2n) 4 He)) = (10) 
1.22 x lO^T-^e- 11 - 327 ^ 1 ' 3 cm 3 s" 1 mole" 1 . 



4. 7 Li( 3 He, 4 He) 6 Li and 7 Be( 3 H, 4 He) 6 Li 

The first reaction was studied at fairly low energies 
[45l ]. Data suggest a fairly constant S-factor for the 
ground state over almost 1 McV in the center of mass. 
Data for the ground state and the first excited state 
(which decays to 2 H+ 4 He) are well established, and an 
approximate total cross section for the second excited 
state (T = 1, so is narrow, and decays to the 6 Li 
ground state) can be inferred. Adding the S-factors for 
7 Li( 3 He, 4 He) 6 Li to the ground state (17.2 MeV-barn) 
and second excited state (12.3 MeV-barn, taken at the 
lowest energy) gives 29.5 MeV-barn. In the absence of 
data for the second reaction, we have assumed that the 
same S-factors for making 6 Li apply to it, as well as for 
destroying 7 Li or 7 Be. The first excited state is excited 
with about 5/7 of the strength of the ground state in the 
7 Li( 3 He, 4 He) reaction. Therefore the assumed rates, all 
in units of cm 3 s _1 mole -1 , are: 

A A (eri;( 7 Li( 3 He, 4 He) 6 Li(ground state))} = (11) 

3.27 x 10 n ^ 2/3 e- M ^ 1/3 , 

iV j4 (cri;( 7 Li( 3 He, 4 He) 6 Li(l st excited states))} = (12) 
1.36 x lO 11 ^ 2 '^- 17 - 98 ^" 3 , 

A A (eri;( 7 Be( 3 H, 4 He) 6 Li(ground state))} = (13) 
= 2.87xl0 n T- 2/3 e- 13 - 722 ^ 1/3 , 

iV j4 (cri;( 7 Be( 3 H, 4 He) 6 Li(l st excited states))} = (14) 
1.19xl0 u T 9 " 2/ V 13 - 722T 9 " 1/3 . 

5. 3 He( 3 H, 7 ) 6 Li 

Despite the difficulty of this experiment (primarily be- 
cause of a huge neutron background), this cross section 
both to the low lying states and to higher lying states has 
been measured [48J . The 6 Li ground state is populated, 
and the second excited state, which also decays to the 
ground state, has a cross section of about 18% of that 
to the ground state. However, the S-factor is difficult to 
ascertain from the data for this reaction. It appears to 
increase linearly from near zero at low energy to a much 
larger value before it levels off [UH^- Since this reaction 
is thought (47T | to result from continuum capture, this is 
surprising. A more likely explanation is that the lower en- 
ergy points distort the results (the beam was stopped in 
the target at the lower energies, complicating the deter- 
mination of the effective target thickness). We therefore 
calculated the S-factor for the points at around 1 MeV 
in the center of mass, and assumed it to be constant, as 
might be expected if the reaction is indeed continuum 
capture. The value for the ground state was increased by 
18% to accommodate the contribution to 6 Li formation 
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from the second excited state. The resulting reaction 
rate, in units of cm 3 s" 1 mole -1 , is: 

A A (a«( 3 He( 3 H, 7 ) 6 Li)) = (15) 
1.02 x 10 7 T- 2/3 e- 7 - 729T o 1/3 . 

6. 7 Be( 3 H,p) 9 Be 

The dominant feature of this reaction appears to be a 
broad resonance in 10 B [4l|, T = 600 keV, at 18.800 MeV, 
seen in several reactions, which is 132 keV above the 
3 H+ 7 Be threshold. Since this is a broad resonance, much 
of which is sub-threshold, it would not be expected to 
contribute strongly to this reaction. Accordingly, in the 
absence of better information, we have simply assumed 
the same rate for this reaction as for the 7 Be( 3 H, 4 He) 6 Li 
reaction (to the excited states). Because this resonance 
would be expected to break up into 4 He+ 4 He+ 2 H most 
of the time (and its decay to 6 Li is included in item 4), 
the rate for this reaction only destroys 7 Be (in units of 
cm 3 s -1 mole" 1 ): 

N A (av( 7 Be( 3 R, a + n) 4 He)) = (16) 
1.19 x 10 n T 9 - 2/ V 13 - 722T 9 " 1/3 . 

7. 7 Be( 3 He,p) 9 B and 7 Be( 3 H,n) 9 B 

These reactions go to a variety of states in 9 B, 
all of which decay by proton emission, including the 
ground state. Since these are deuteron transfer reac- 
tions, and have comparable (large) Q-values to that 
for the 7 Li( 3 H,n) reaction, we assumed the same astro- 
physical S-factor for it to the 9 B ground state as for 
7 Li( 3 H,n) 9 Be(ground state) [H (including more than the 
leading term does not seem justifiable), and noted 
that the 7 Li( 3 He,p) 9 Be yield to all states appeared to be 
greater than that to the ground state by approximately 
a factor of five. After appropriate corrections for the nu- 
clear charges [38}, the assumed rates are 

N A {av( 7 Be{ 3 Be,p) 9 B)) = (17) 
2.06 x 10 u T 9 - 2/3 e- 21 - ra2T »~ 1/3 , 

N A (av{ 7 Be( 3 R,n) 9 B)) = (18) 
1.64 x 10 n T 9 - 2/ V 13 722T 9 ~ V3 , 

where N A (av) is in the usual units of cm 3 s -1 mole -1 . 
Another reaction, 9 Bc( 3 Hc, 6 Li) 6 Li, was not included in 
the reaction network, as it was found 54] to have a much 
smaller cross section than 9 Be(p, 4 He) 6 Li, and it depends 
on the 9 Be abundance, which is small in BBN. Thus it 
would not be expected to contribute appreciably to 6 Li 
abundance, despite the fact that it makes two 6 Li nuclei, 
or to 9 Be destruction. 



8. 7 Li( 3 He,d) s Be, 7 Li(d,n) 8 Be, 7 Be(d,p) 8 Be, and 
7 Be( 3 H,d) 8 Be 

These reactions are all single nucleon transfers, and all 
have large, and comparable, Q-values. We therefore esti- 
mated an S-factor for the 7 Li(cZ, n) 8 Be reaction from the 
total cross section for that reaction [43[ , and applied it to 
all four reactions. That S-factor is reasonably well rep- 
resented by a constant of about 2.5 xlO 4 keV- barns (and 
which is lower than, but reasonably consistent with, that 
from Angulo et al. [Ill, measured for 7 Be(d, p)2 4 He); we 
therefore assumed that value, and then calculated the 
reaction rate from the standard expression |38| . includ- 
ing also spin statistics. Note that the 7 Li(<i, n) 8 Be reac- 
tion proceeds almost entirely through two states at 16.63 
and 16.92 MeV in 8 Be for energies in the compound sys- 
tem more than 1 MeV above the energy of those states 
55]. However, the Q-values for all except 7 Be(d,p) 8 Be 
are well below that energy, and even its Q-value is just at 
that suggested threshold [56|. Thus the assumed rates, 
in units of cm 3 s -1 mole -1 , are: 

N A {av( 7 U( 3 lle,e) s Be)) = (19) 
6.26 x 10 11 T 9 - 2/3 e- 17 - 980T »" l/> , 

N A {av( 7 Li(d,n) s Bc)} = (20) 
2.43 x lO 11 ^ 27 ^ -10 - 254 ^ 173 , 

N A (av( 7 Be{d, P ) 8 Be)) = (21) 
2.67 x 10 11 T 9 " 2/3 e- 12 - 422T <T 1/3 , 

N A (av( 7 Be( 3 R,d) 8 Be)) = (22) 
5.47 x lO 11 ^" 27 ^ -13 - 722 ^ 173 . 

A recent paper by Cyburt and Pospelov [32j suggests 
that there might be a resonance in the 7 Be+ 2 H reac- 
tion that would decay to a variety of final states, all of 
which would destroy 7 Be, and that this might solve the 
problem of strong 7 Be production in BBN. This asser- 
tion was based on a paper by Dixit et al. [36], which 
studied 9 Be(p,p), and which identified a strongly ex- 
cited, and relatively narrow, state at 16.7 MeV. In that 
paper the possible structure of a narrow state at such 
high excitation is discussed, and is identified with earlier 
theoretical work [37] in which this state is suggested as 
having a structure of an Si/ 2 neutron coupled to highly 
excited 2 + core states in 8 Be. Such states could only 
be v ld v p 2 -> n sd 7r p i ! ) or VsditsdVp 1 tt~ 1 configurations, so 
would require a complicated (and therefore weak) multi- 
step reaction to excite any reactions involving 2 H+ 7 Be, 
since 7 Be presumably has a dominant configuration of 
Up on an 8 Be core. As a consequence, it seems unlikely 
that reactions proceeding through this state could solve 
the 7 Be excess problem. Nonetheless, the considerable 
leverage that this reaction can have on 7 Be, as pointed 
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out by Cyburt and Pospelov |32| . remains alluring. 



9. 7 Li( 3 He,n) 9 B 

Data for this reaction are sparse [42[ , but suggest that 
many states in 9 B are populated. A spectrum at 2 MeV 
suggests that the strength to the ground state of 9 B con- 
stitutes roughly 20% of the yield (which, because the 
Q-value is so large, we took to be constant with energy). 
A rough S-factor can be obtained from the data for the 
transition to the ground state; we found it to be 180 MeV- 
barn, although this determination is complicated by an 
apparent resonance structure above the region of inter- 
est and a lack of information about experimental details. 
This gives a reaction rate, in units of cm 3 s _1 mole -1 , of 



N A {crvCU{ 3 Re,e) 8 Be)) = 
2.0xl0 12 T 9 2/3 e- 17 - 980T 9 1/3 . 



(23) 



1e-05 



1e-25 



N/H 




10. 7 Be(d, 3 He) 6 Li and 7 Be*(d, 3 He) 6 Li 

There are no data for the first reaction 7 Be(e?, 3 He) 6 Li. 
The Q-value for this reaction is slightly negative. How- 
ever, if 7 Be is in its first excited state at 0.429 MeV, 
which would be populated to some extent in the high 
temperature thermal environment in which BBN occurs, 
this situation changes. Then a strong resonance, with 
a width of 22 keV, would be expected at an excitation 
energy of 17.190 MeV in the 9 B compound nuclear sys- 
tem, which is seen in the 6 Li( 3 He,p) reaction [57[. In the 
absence of information about the partial widths of this 
state, we have assumed, in order to obtain an upper limit 
on its effect, that the partial widths for the 6 Li+ 3 He and 
7 Be+ci channels were equal at 11 keV. Since the spin of 
the resonance is also not known, we have taken the spin 
statistical factor in Eq. (Q| to be 1.0. Then the reaction 
rate for this reaction is found to be: 



AUW 7 Be*(d, 3 He) 6 Li)) = 
4.31 x 10 8 T- 3/ V 18336 / T 9, 



(24) 



where Na(<jv) is in the usual units of cm 3 s -1 mole -1 . 
In the BBN reaction network and abundance calculations 
the first excited state 7 Be* was entered as a separate 
nuclear species. The abundance of this species was tied 
to the abundance of the ground state by multiplying by 
the thermal population factor 

_ (2 J* + 1) ■ e- E -' T 2 • e -4-98/T 9 

"* — v ~ A ' (■ 5 ' 

^Be 4 

where J* = 1/2 and E* w 0.4292 MeV are the 7 Be first 
excited state spin and excitation energy respectively, and 
where we approximate the 7 Be nuclear partition func- 
tion as Ztq c r< 2J groun d + 1 = 4. At Tg « 1, where 
7 Be is abundant enough to be affected appreciably by 
nuclear destruction reactions, this population factor is 
small, P,k3x 10~ 3 . 



FIG. 2: Nuclear abundance as a function of temperature Tg, 
where Tg = T/10 9 K. Abundances are given as mass fraction 
for 4 He and number abundance relative to hydrogen for all 
others. 



III. BBN RESULTS FROM ADDED NUCLEAR 
REACTIONS 



We have employed a variant of the standard BBN re- 
action code @. Our code is described in detail in Refs. 
pll [l2| • The reaction network, which includes the stan- 
dard model nuclei and reactions plus the reactions we 
have added, is indicated in Fig. [TJ Note, however, that 
this cartoon does not show all the reactions in the net- 
work, but rather is simply a guide to some of the new 
reactions included here, including the inclusion of 6 He to 
the network. The reactions that were added to the BBN 
code, along with their Q-values, are listed in Table I. 

The evolution of the abundances in our calculations is 
shown in Fig. 2. The importance of the various reactions 
to the BBN abundances are specified below and indicated 
in Table 1. In general, it can be concluded that the added 
reactions have very little effect on BBN. However, since 
some of the reaction rates were quite uncertain, we in- 
creased the more uncertain rates by a factor of 1000 to see 
if such huge increases in the rates might have an effect. 
In most cases, only tiny effects were observed; the des- 
ignation "none" means that the BBN abundances were 
not changed by more than 0.1% with the factor of 1000 
enhancement for that rate. In this table, "small" means 
there was an alteration in abundance yield by more than 
0.1% when the rate was enhanced by a factor of 1000. 
Likewise, in this table "large" means a 10% or bigger 
abundance change when multiplying our best guess rate 
by 1000. 

It is remarkable that virtually nothing can be done in 
the context of the nuclear reactions to increase the 6 Li 
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FIG. 3: 6 Li destruction rates as a function of temperature Tg. 



abundance up to the value suggested by the observations, 
even with the huge multiplicative factors used, although 
the 3 H( 3 He,7) 6 Li reaction might warrant further exper- 
imental study. Only the 7 Bc( 2 H, 3 He) 6 Li reaction pro- 
duced a large enough effect to suggest that it might be 
important in destroying some of the mass-7 nuclei ul- 
timately resulting from BBN. However, that required a 
factor of xlOOO increase in the reaction rate, so that effect 
must be viewed with skepticism. 

Of some interest was the inclusion of 6 He in the re- 
action network. This species is converted quickly to 6 Li 
by weak interactions. Although the 6 He total weak de- 
cay rate is accelerated by v & and e + capture processes at 
high temperature, by Tg ~ 1 the laboratory beta decay 
rate (~ 0.9 s _1 ) dominates [58|. Despite the rapid weak 
decay, as can be seen in Fig. [51 a significant 6 He abun- 
dance is built up during BBN. Although this 6 He produc- 
tion initially contributes to 6 Li through decay and lepton 
capture reactions [58(, the 6 Li destruction processes ef- 
fectively demolish any potential increase in the final 6 Li 
abundance. All of the 6 Li destruction reactions included 
in the BBN code are shown in Fig. [3] Of these processes, 
we see that 6 Li(p, a) 3 He has the largest effect in 6 Li de- 
struction mainly because of the abundance of protons 
when 6 Li is created. 

Following the Cyburt-Pospelov [32j suggestion on the 
sensitivity of the 7 Be destruction to new nuclear physics, 
we performed simulations of BBN in which the rate for 
7 Be(d,p) 8 Be was used as a surrogate for all the deuteron- 
induced reactions that destroy 7 Be. We increased that 
rate by factors of 10 and 100 from the standard rate. The 
factor of 10 produced little change in the 7 Be yield, while 
the factor of 100 did reduce the yield by 30%, but still 
not enough to achieve the required reduction factor of 
roughly 3. From our investigations, we believe that even 



the factor of 10 increase in this rate is unrealistically 
large. 



IV. NON-THERMAL NEUTRONS IN BBN 

In order to proceed with the analysis involving the 
non-thermal neutrons, it is essential to describe their en- 
ergy spectrum. The 14.1-MeV neutrons produced by the 
3 H((i, n) 4 He reaction in BBN lose energy and eventually 
thermalize by scattering with other nuclei. If it is as- 
sumed that the time in the big bang is late enough that 
scattering from other neutrons can be ignored, the most 
important scattering target is -"-H, which makes up over 
90% of the number density and is also the lightest iso- 
tope. The only other relevant nuclide, 4 He, has a neutron 
scattering cross section w 30% larger than 1 H, but con- 
tributes less than 10% of the total number density and 
is also less favorable for energy transfer due to its larger 
mass. Thus we only consider energy loss due to scattering 
from 1 H. Our ultimate goal is to calculate the reactions 
induced by these neutrons before they thermalize. 

Data for the n — p total cross section are given in 
Ref . [59T - [6l| ; we utilize these along with a power-law pa- 
rameterization 



a(E)=a {E/E ) a , 



(26) 



where E is the neutron energy and the constants are 
given by <tq = 0.685 barns, Eq — 14.1MeV, and a = 
—0.834. This parameterization reproduces the experi- 
mental data within 10% for 3 < E < 20 MeV. It over- 
estimates the cross section for lower energies. However, 
we will use this parameterization in order to produce an 
upper limit on the effect of the non-thermal neutrons. 

The energy transfer to the proton depends on the neu- 
tron scattering angle. Assuming that (1) the n-p differ- 
ential cross section is isotropic in the center-of-mass sys- 
tem, (2) the neutron mass is equal to the proton mass, 
and (3) non-relativistic kinematics are valid, the energy 
loss of the neutron due to a single scattering event is uni- 
formly distributed, from zero up to the initial neutron 
energy. These conditions are satisfied within a few per- 
cent for neutrons below 20 MeV; this assumption leads 
to significant simplifications in the analysis and will be 
utilized below. 

We assumed that neutron thermalization occurs on a 
time scale much shorter than the local Hubble time and 
the timescale for changes in nuclear abundances. Then 
the probability for neutrons of energy E scattering from 
hydrogen in a time interval At is given by 



P{E) = <7(E)n u vAt, 



(27) 



where cr{E) is the n-p total cross section, nn is the : H 
number density, v is the neutron speed, and At has been 
assumed to be sufficiently small that P(E) « 1. Like- 
wise, the probability for neutrons of energy E to scatter 
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from hydrogen into the energy interval [E',E' + dE'] is 
given by 



P(E,E')dE' = criE^nvAtdE', 



(28) 



where we have made use of the uniform energy distribu- 
tion discussed above. 

The non-thermal neutron energy distribution consists 
of two contributions: a mono-energetic component of no 
neutrons per unit volume accounting for 14.1-MeV neu- 
trons which have not yet scattered, and a continuous dis- 
tribution f(E) describing the number of neutrons per 
unit volume per unit energy at lower energies. We define 
E = 14.1 MeV. 

f(E) will evolve over the time interval At as: 

A/(£) = -f(E)P(E) + n P(E ,E) (29) 



f(E')P(E', E)dE' . 



The first term on the right hand side of Eq. (1291) de- 
scribes the loss due to downscattering to lower energies, 
the second describes the increase due to the downscat- 
tering of 14.1-MeV neutrons, and the third describes the 
increase due to the downscattering of continuous neu- 
trons at higher energies. With the assumption of a fast 
timescale for thermalization, f(E) will be in dynamic 
equilibrium and Af(E) — 0. Using these equations and 
E = (l/2)mv 2 gives 



1/2 



f(E)a(E)E L ^ +n a(E)E 
f(E , )o{E'){E')- 1 l 2 dE' = 0. 



(30) 



This implies that 



f{E ) = n /E . 



(31) 



If <j(E) is given by Eq. (f2l)|). Eq. (f3"Tj) can be solved ana- 
lytically to yield the solution (for a ^ —1/2), 



f(E) 



no 



E(a + 1/2) 



{E /E) a+1 / 2 + a- 1/2 . (32) 



This distribution for f(E) is shown in Fig. |4] Because 
our model does not include any neutron sinks or upscat- 
tering, the distribution is highly singular at E = 0. The 
normalization down to some minimum energy can be de- 
fined: 



N m =^lf{E)dE=j^ JW (33) 
x [(E /E m r +1 / 2 - 1 + (a 2 - 1/4) log(E /E m )} . 



A minimum energy of E m —2.5 MeV gives N rn — 2.975no. 

The number density of 14.1 MeV neutrons no can be 
found by equating the rates for production and destruc- 
tion. The destruction rate is given by: 




FIG. 4: Neutron energy distribution. 

where Eq = (l/2)mv 2 . The production rate is given by 

r = n t n d {(Jv)td^na (35) 

and consequently no can be calculated using 

n = [n t n d {av) td ^ na ]/[n H a(E )v }. (36) 

Note that the quantities njj, n d , n t and (crv)t d -+na are 
available in existing standard BBN calculations. 

The reaction rate (reactions per unit volume per unit 
time) of the non-thermal neutrons with species x is given 
by 

r x = n x n a- x (E )v +n x / f(E)a x (E)vdE, (37) 

where a x [E) is the n + X cross section. In this model 
the only temperature dependence is through n . A 
temperature-independent reaction rate {a x v} can thus be 
defined as 

(a x v) = {a x {E )v ) + (1/n) [ ° f(E)a x (E)vdE (38) 

J E m 

which only needs to be calculated once. The rate of reac- 
tions defined by Eq. (|34[) can now be written 



r x = n x n {a x v). 



(39) 



r = n P(E) = n Q n H ((j(E )vo}, 



(34) 



With the above framework, it is now simple to implement 
the extra reaction flow due to non-thermal neutrons into 
the BBN calculations. 

These show that the neutron energy distribution for 
high energy neutrons, that is, those above 2.5 MeV, is 
many orders of magnitude larger than the prediction 
of the Maxwell-Boltzmann distribution when the non- 
thermal contribution is included. See figures 0] and [S] 
The effects on the BBN abundances, however, were found 
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FIG. 5: The fraction of neutrons with energies above E m = 
2.5 MeV is shown by the solid curve. The fraction predicted 
by Maxwell-Boltzmann distribution is shown by the dashed 
curve. 




20 



neutron energy (MeV) 



FIG. 6: The cross section for the n + d — > 2n + p reaction as 
a function of incident neutron energy, from the ENDF/B-VI 
evaluation [sil ] . 

to be at the 10~ 4 level or less, much smaller than the level 
of uncertainties resulting from other sources such as reac- 
tion rates. Details of several of the reactions considered 
are given below. 

11. 2 ¥L(n,p)n 

Deuterium formation is crucial to all subsequent nu- 
cleosynthesis. Non-thermal neutrons might allow the en- 
dothermic deuterium destruction 2 H(n,p)n reaction to 
occur at an appreciable rate once nuclear statistical equi- 



librium for n + p ^ 2 H + 7 breaks down following alpha 
particle formation. The cross section for 2 K(n,p)n is 
shown in Fig. |51 The inclusion of this reaction is found 
to increase the deuterium abundance by only 0.011%, and 
the effects on other isotopes are an order of magnitude 
smaller. Note that this is an upper limit on the possible 
effects of the non-thermal neutrons on this reaction. 



12. 3 He(n,p) 3 H 

This reaction is exothermic, so can also proceed also 
via thermal neutrons on this reaction. Thermal neutrons 
are far more effective than non-thermal neutrons in pro- 
ducing reactions as a result of the 1/v law dependence 
of the neutron capture cross section. We therefore con- 
clude that non-thermal effects will be negligible for this 
reaction. 



13. 7 Be(n,X) reactions 

Neutron-induced reactions on 7 Be producing the fol- 
lowing exit channels were considered: 2 H+ 6 Li, n+p+ 6 Li, 
n+ 3 He+ 4 He, and 4 He+ 4 He. The first three are en- 
dothermic; the latter is exothermic, but highly sup- 
pressed for low-energy neutrons, as it cannot proceed 
via s-waves in the entrance channel. The cross sections 
for the reactions are not known. The maximum inelastic 
cross section for n+ 7 Li approaches 3 barns (62|. We have 
adopted this as an upper limit for the above reactions. 
Inclusion of 7 Be(n, n 3 He) 4 He, then, was found to reduce 
the 7 Be abundance by 0.015%. Results for the other exit 
channels would be expected to be similar. Again, these 
are expected to be upper limits on the effect of the non- 
thermal neutrons. 



V. FURTHER DISCUSSION OF THERMAL 
POPULATION OF EXCITED STATES AND 
RESONANT ENHANCEMENT 

The lowest-lying first excited states of the lightest nu- 
clei are those of the mirrors 7 Li and 7 Be, at excitation 
energy 478 keV and 429 keV, respectively. These en- 
ergies are tantalizingly close to BBN thermal energies, 
for example with alpha particle assembly at kT ^100 
keV. Accordingly, we have examined how potentially 
key nuclear reactions like 7 Be(d, 7) 9 B, 7 Be(c£, p)2a, and 
7 Be(<i, 3 He) 6 Li could be altered if they were to proceed 
through a thermally populated first excited state, 7 Be*. 
We have discussed this at length above (section 11.10) for 
7 Be(<i, 3 He) 6 Li. Some further remarks on possible ther- 
mal excitation effects are in order. 

Since the compound nucleus for 2 H+ Be, 9 B, has a 
ground state that particle decays, and this reaction is 16 
MeV up in excitation in the compound nucleus, it is un- 
likely that it would be strongly resonant. Without the 
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sort of enhanced rate that would be produced by a res- 
onance, it could only have a very small effect on BBN 
abundances. Thus, we did not include the 7 Be( 2 H,p) re- 
action in the BBN code. The 7 Be*(p, 7) 8 B reaction ini- 
tiated from the thermally excited 7 Be first excited state 
could produce an enhanced rate. However, it cannot have 
much of an effect on the BBN abundances. There is a 40 
keV wide state in 8 B, the compound nucleus, at 770 keV 
excitation, which from the 7 Be first excited state, would 
be at 241 keV bombarding energy. However, the spins 
and parities of the entrance channel particles together 
with the spins and parities of either the ground state of 
8 B or this resonant state (2 + and 1 + respectively) require 
a p-wave capture. Furthermore, this is a radiative cap- 
ture, for which the cross section is generally roughly two 
orders of magnitude suppressed from that of a particle 
transfer reaction. When one imposes a Boltzmann fac- 
tor to represent the abundance of the excited state, this 
reaction is further suppressed. Consequently we did not 
include this reaction in the BBN code. 

In the end, however, standard BBN yield alterations 
due to reactions through the 7 Be* state are found to be 
small on account of its small thermal population. 

VI. CONCLUSION 

We have examined a number of nuclear reactions likely 
to occur in BBN. Some of these have not been included in 
BBN calculations before and some might be thought to 
figure in the 6:7 Li problems. We have discussed nuclear 
physics uncertainties and other issues (e.g., non-thermal 
neutrons, thermal excitation of nuclear target states) and 
we have performed sensitivity studies with an extended 
BBN reaction network calculation of abundance yields 
to get an idea of what might be possible with these reac- 
tions. 

In broad brush, the reactions added to the BBN code 
had virtually no effect on the BBN abundances, although 
a few of them did have an effect when their rates were 
increased by the factor of xlOOO. In particular, the 
3 He( 3 H,7) 6 Li reaction did produce a large 6 Li abundance 
enhancement, and the 7 Bc( 2 H, 3 He) 6 Li reaction produced 
a 30% decrease in the mass-7 abundance when their cor- 
responding rates were increased by a factor of 1000. Al- 
though this multiplier is probably well outside the ex- 
pected uncertainties for the rates for these reactions, 
these reactions might bear further study just to pin down 
the BBN abundances with the highest possible precision. 

More important, however, is the conclusion from this 
study that the solutions to the problems of the 7 Li and 
6 Li abundances are hard to find in the reaction rates. 
For 6 Li, it may well be that the observed BBN abun- 
dances are simply incorrect; it has been suggested by 
Cayrel et al. [27J that the claimed 6 Li abundance is ac- 



tually a spurious effect that is generated by convective 
Doppler shifts. However, the Li discrepancy is not as 
easy to ignore. If both problems persist, however, their 
solution may well require nonstandard physics. Several 
papers have been published recently, some involving un- 
stable particles in the early universe that could have af- 
fected BBN. One set of papers studied the effects of the 
reactions resulting from non-thermal decay products of 
the particles. These could destroy 7 Li directly, reduc- 
ing its abundance, and would make 6 Li by producing 
non-thermal distributions of 3 H and 3 He from proton 
and neutron knockout reactions on 4 He. These could 
then undergo 4 He( 3 H,n) 6 Li and 4 He( 3 He,p) 6 Li reactions 
respectively. However, as discussed by Kusakabe et al. 
[631 ] . these solutions have the difficulty of also affecting 
the 2 H and or 3 He abundances, creating new problems. 
Note, however, that new cross sections for 4 He(7, p) 3 H 
and 4 He(7,n) 3 He may affect that conclusion J64|. The 
second possible effect of unstable particles in the early 
universe would be if they were massive and negatively 
charged, so could attach to the nuclei as they were pro- 
duced. A suggestion by Pospelov '65], that such particles 
could enhance the 4 He(<i, 7) 6 Li reaction by many orders 
of magnitude by enabling a transfer reaction instead of 
the radiative capture is the critical feature of this solu- 
tion. This has been found by several authors to solve 
both Li problems [MlMS]. 

In summary, we find little effect on standard BBN 
abundance yields with the addition of the new nuclear 
physics we discuss here, even given a fair uncertainty in 
issues that bear on key reaction rates. It could be argued, 
however, that the new nuclear reaction physics we exam- 
ine here may yet be important for non-standard BBN 
scenarios with new particle physics. Indeed, a conclusion 
of this paper is that there seems to be little chance of 
solving either of the "lithium problems" by conventional 
nuclear physics means and, if these problems stand up 
to future observations, we may be forced into just such 
non-standard BBN scenarios. 
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